Thermoelectric performance of a driven double quantum-dot 
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In this paper we investigate the thermoelectric performance of a double-dot device driven by time- 
dependently modulated gate voltages. We show that if the modulation frequency f2 is sufficiently 
small, not only quantized charge pumping can be realized, but also the heat current flowing in the 
leads is quantized and exhibits plateaux in units of kB.T\n2. The factor ln2 stems from the 
degeneracy of the double-dot states involved into transport. This opens the possibility of using the 
pumping cycle to transfer heat against a temperature gradient or to extract work from a hot reservoir 
with Carnot efficiency. However, the performance of a realistic device is limited by dissipative effects 
due to leakage currents and finite-frequency operation, which we take into account rigorously by 
means of a real-time diagrammatic approach in the regime where the double dot is weakly coupled to 
the leads. We show that despite these effects, the efficiency of a double-dot charge pump performing 
work against a dc-source can reach of up to 70% of the ideal value. 

PACS numbers: 73.63.-b,73.23.Hk,73.50.Lw,85.80.Fi 



I. INTRODUCTION 

The tunability of nanoscalc systems such as quantum 
dots or single-electron boxes allows to exploit their func- 
tionalities in the realization of nanoscalc electric and 
thermoelectric devices. While the use of quantum dots 
and metallic islands as single-electron transistors^ ultra- 
sensitive detectors^ or thermometers^ has already been 
considered for a long time, it is only more recently 
that their applications as thermoelectric engines^— 
refrigerators and heat rectifiers^— - to name some 
examples - started to attract considerable attention. 
Since Coulomb interaction and quantum interference ef- 
fects play an important role in these devices, their ther- 
moelectric properties differ strongly from those of clas- 
sical macroscopic systems ^Sr— and can even result in a 
significant increase of the thermoelectric efficiencyi 2 *^— 
Moreover, thanks to their stability and to the possi- 
bility of monitoring their state with integrated charge 
or temperature detectors, quantum dots and single- 
electron boxes represent an ideal playground to experi- 
mentally test the predictions of so-called fluctuation the- 
orems; 2 ^— or to investigate possible implementations of 
Maxwell's demons in solid state environment i 32 1 33 

Our particular interest concerns the operation of a 
quantum dot system driven by time-dependent fields. 
Since the original experiment of Pothier et al£^ it is 
well known that by changing slowly the gate voltages 
applied to a double-dot (or a double metallic island as 
used in Ref. H3), electrons can be transferred from one 
lead to the other in a controlled way, generating a quan- 
tized dc-currcnt / = ±ef2/27T, where Q is the frequency 
of the ac-signal applied to the gates. This result has 
a clear interest for metrology^ and stimulated an in- 
tense experimental 3 ^— and theoretical activit y 40 ' 41 Im- 
portantly, single-electron pumping can be achieved both 



when the two leads are in the same equilibrium state and 
when a finite bias voltage is applied to the device, as well 
as it is possible to pump electrons against the preference 
direction set by the bias. 

Beyond the goal of realizing quantized charge pump- 
ing, there has been interest in time-dependently driven 
double-dot systems thanks to their complex internal 
structure, which results in highly sensitive devices. To 
name some examples, time-dependently driven hybrid 
double-dot devices with ferromagnetic leads have been 
shown to allow for pure spin currents,— ~— the orbital 
degree of freedom of a double-dot allows to study spin- 
orbit coupling effects in pumping ) 45 ' 46 and a double-dot 
pump with superconducting leads has been proposed as 
a detector for crossed Andreev reflection.— 

Motivated by the possibility of pumping charge against 
a bias, we study the performance of a driven double- 
dot as a nanoscale "battery charger" transferring elec- 
tron from a lower to a higher chemical potential, as 
well as we investigate its efficiency as a heat pump (or 
heat engine) operating between two reservoirs at differ- 
ent temperatures. To this end, we consider the charge 
and heat currents flowing in the leads in response to 
the time-dependent driving. Heat currents in electronic 
quantum-pumps have been so far mostly investigated 
in the limit where the electronic interactions can be 
neglected^ ' 10 ' 48- — The operation of the double-dot pump 
is however based on charging effects^ and therefore we 
adopt a real-time diagrammatic approach^ which al- 
lows to take into account arbitrarily strong Coulomb in- 
teraction between electrons in the dots, as well as non- 
equilibrium conditions induced by dc voltage and temper- 
ature gradients, and the time-dependent driving] 53 ' 54 In 
the limit of weak coupling to the leads and slow driving, 
we perform a rigorous expansion in the modulation fre- 
quency Q, taking into account contributions to the charge 
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FIG. 1. (Color online) (a) Sketch of the potential landscape of 
a double-dot setup, where time-dependent gates modulate the 
positions of the dot levels, (b) Stability diagram of the double- 
dot for the case V — 0: black full lines indicate the borders of 
the stability regions for negligible inter-dot coupling, while the 
red dashed lines show the resonance positions of bonding and 
antibonding state, assuming U — 20£;bT and t c = lOfceT. 
The orientation of the system of coordinates in parameter 
space formed by the detuning e — el — £r and the mean 
energy E — (el + £r)/2 is also shown, (c) Sketch of possible 
pumping trajectories in parameter space. The color-scale plot 
represents the dc charge-current through the double-dot in 
the linear response regime (eV = 0.5fcBT); (d) Dc charge- 
current at finite bias (eV = lOksT). Other parameters in 
panels (c) and (d): U = 20fc B T, U = l0k B T, T L = T R = T/2, 
hT = fc B r/4 and TL = Tr = T. 



and heat currents in zeroth order in the frequency (dc 
leakage-currents), in first order in Q, (genuine adiabatic 
pumping), as well as in second order in the frequency, 
which arc particular important in order to account prop- 
erly for the heating effects due to the ac- driving. 

We show that in certain regimes not only the pumped 
charge-current is quantized, but also the heat current 
shows well defined plateaux which are directly related 
to the electronic temperature in the leads and the degen- 
eracy of the double-dot states involved into transport. 
Moreover, we show that in these regimes the pumping cy- 
cle can be considered as a close analog of the Carnot cycle 
and that the double-dot pump can in principle be em- 
ployed to transfer heat against a temperature gradient, 
to extract work from a hot reservoir, or to move charges 
from a lower to a higher chemical potential with maxi- 
mal efficiency, if the driving is infinitely slow. The per- 
formance of a realistic pump is however limited by leak- 
age currents and heat production due to finite-frequency 
operation. We investigate these limitations in detail and 
find that, at least for the case of the charge pump working 
against a dc-bias, efficiencies up to 70% of the maximal 
value can be obtained. 



The manuscript is structured as follows. We intro- 
duce the model for the double-dot and the real-time- 
diagrammatic technique used for the calculation of the 
charge and heat currents in Section [TTJ In Section IIIII we 
discuss first the case of pure adiabatic pumping of charge 
and heat, and then present the operation of three ideal 
double-dot based engines. Finally, the limitations to the 
perfect operation of these engines due to leakage currents 
and finite- frequency driving are discussed in Section [IVI 



II. MODEL AND TECHNIQUE 
A. Double dot in the molecular regime 

We consider a double-dot device formed by two single- 
level spin-degenerate quantum dots connected in series 
and coupled to external leads, see Fig. [TJ,. This system 
is described by the Hamiltonian 



H — Hdd + E H a + H t 

a— L,R 



(la) 



Here, H^d is the Hamiltonian of the isolated double-dot: 

TTl 

Hdd = e 

I ~ a— L.R 

(lb) 



= ^ + UnhfiR + — ^ n a (h a - 1) 

a— L,R a— L,R 

tc 

2 



<7=t4 

The dot operators 4,r(4(r) create (annihilate) an elec- 
tron with spin a =f, 4- an d energy e a in the dot a = L, R. 
The corresponding number operator of electrons in each 
dot is given by n a = h at + h ai with h aa = cP aa d aa , 
we also define the total occupation-number operator of 
the double-dot, h — hh + riR. The inter- and intra-dot 
Coulomb interaction is denoted by U and U' , respec- 
tively. In the following, we will assume the onsite in- 
teraction U' to be the largest energy scale in the sys- 
tem (U' — > oo), so that each dot can be at most singly 
occupied. Hopping from one dot to the other occurs 
with the inter-dot coupling amplitude —t c /2, where t c 
is taken to be real and positive. The single-particle ener- 
gies e a = e a (t) of the dots can be tuned by external gates 
locally applied to the two dots and are in general time- 
dependent. The leads are described as non-interacting 
Fermi liquids with the Hamiltonian 



H a = 



e akaC i aka C aka , 

fc,(T=t,| 



(lc) 



at 



r (c Qfccr ) are the creation (annihilation) opera- 
tors for an electron with momentum k and spin a in lead 
a. Finally, the coupling between the double dot and the 
leads is given by 



Htun = / ] 



(Id) 
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The coupling is quantitatively characterized by the 
energy independent tunnel-coupling strength r a = 
2irv a \t a \ 2 /h. Here t a is the tunneling amplitude between 
dot a and its neighboring lead, which we assume to be 
spin and momentum independent, and v a is the density 
of states of lead a in the wide band limit. 

The imposed boundary conditions are such that each 
lead is in local equilibrium with the electrochemical po- 
tential [i a and the temperature T a . Here we take the 
mean electrochemical potential /i = ^ L + Mfl as the zero 
energy level, i.e. /i — 0. All single-particle energies in 
Eq.JTJ) are expressed with respect to it. The difference 
between the electrochemical potentials of the two leads 
is fixed by the applied bias voltage V. In the rest of the 
paper we will assume /il = — fJ-n = eV/2 > 0, with e > 
the absolute value of the electron charge. 

In the following we will focus on the case of strong 
inter-dot coupling, namely when the inter-dot dynamics 
is much faster than the dot-lead hopping (t c 3> HT with 
r = Tl + Tr). In this limit, it is useful to diagonalize the 
single-particle sector of Hdd introducing its "molecular" 
eigenstates, namely the bonding and antibonding states 
| be) = d^lO) and |aer) = d a(T |0), which go along with the 
creation operators 



are defined by the conditions = and e a + U = and 
an anti-crossing behavior is shown at the triple points. 
This anti-crossing can be clearly distinguished in the dc- 
current in the linear response regime (see Fig. [TJ;), and 
is a hallmark of a tunnel-coupled double-dot system. For 
finite bias the triple points get broadened (see Fig. [T}1) , 
as there are now entire regions in parameter space that 
correspond to no stable charge configuration. 



B. Kinetic equation 

We want to describe the dynamics of the double dot 
in the presence of a time-dependent driving applied to 
its gates in the limit of weak coupling to the leads 
fir <C ksT. The state of the double-dot at a given 
time t is described in general by its reduced density ma- 
trix pdd{t)- However, if the inter-dot coupling is much 
stronger than the one to the leads (t c ^> HT), the dy- 
namics of the diagonal and off-diagonal elements of pdd it) 
decouple in lowest order in T, and we can restrict our- 
selves to study only the diagonal elements, namely the 
occupation probabilities.— Their evolution is governed 
by the generalized master equation^ 



dt, = 

b/a a 



V2 



'1±- 



'it 



*Lcr 



These states have the energies 



1 



e b/a 



(2) 



where the upper (lower) sign corresponds to the bonding 
(antibonding) state and E = (cl + £r)/2 and e = £l — 
£r are the mean energy and the detuning between the 
energies of the two dots, respectively (see Fig. [Tp). In 
this new basis, the tunnel coupling between the double- 
dot and the two leads is expressed by effective rates for 
tunneling via the hybrid states |ber) and |atr)4^ 



(3a) 
(3b) 

We also define the total broadening of the bonding and 




antibonding states as T 



b/a 



/a 



Because 



of the dependence on the level detuning e, these effective 
rates become a function of time whenever an out-of-phase 
modulation is applied to the gates of the two dots. 

The stability diagram of the double dot is shown in 
Fig. [Tp, and identifies regions of different equilibrium oc- 
cupation numbers for the two dots, as a function of the 
energies cl and £r. The points where three charge states 
are degenerate are named triple points. For strong inter- 
dot coupling t c ^> HT, the edges of the stability regions 



dt 



Pit) 



dt'W(t,t')p(t'), 



(4) 



with p = (po , Pat , Pai , Pbt , Pb4. , Ptt i Pti > Pit ' Pu) i where 
we omitted the time arguments for simplicity. These are 
the probabilities that the double dot is empty, po, that 
the double dot is singly occupied with an electron with 
spin f or J, in the bonding or in the antibonding state pbo- 
andp a(T , or finally that the double-dot is doubly occupied, 
p CTCT ', where the two subscripts indicate the spin of the 
electrons occupying the left and right dot, respectively. 
As mentioned before, double occupation of a single dot 
is energetically forbidden due to strong onsite repulsion. 
The kernel W(£, t') is a transition matrix that describes 
tunneling events due to the coupling to the leads. It has a 
functional dependence on the time-dependent parameters 
e a (t). 

We focus here on the case of slow driving and weak 
coupling to the leads. In this regime, Eq.(|3]) can be solved 
perturbatively by performing a double expansion in the 
driving frequency ft and in the tunnel coupling This 
gives rise to a hierarchy of equations that in the limit of 
weak coupling to the leads {HT fcgT) reduces to 53 ' 54 



WtPl 



(o) 



0. 



WtP { t k) - Jt P\ 



(fc-1) 



(5a) 
(5b) 



where the superscript (k) indicates terms of order fi fe 

in the expansion p = X)fc>oP*^ • ^ ne ma trix Wt is 
the zero-frequency Laplace transform of the Kernel in 
time-space, Wt = /_ d(t' — t)~W(t,t'), and it contains 
the transition rates between different double-dot states 



4 



to first order in V (i.e. as given by Fermi golden rule), 
evaluated with all time-dependent parameters "frozen" 
at time t. It satisfies the sum rule X^JWtlu = 0, which 
ensures the conservation of probability, and can be writ- 
ten as a sum of independent contributions from the two 
reservoirs W t = S a =i,R W «,*- 

Together with the normalization conditions 

J2iPi°t = 1 an d J^iPi^ ^ = 0' Eqs.JSJ allow to evaluate 
iteratively all p t , given the knowledge of the instanta- 
neous solution p[°'££ The latter represents the steady 
state the system would relax into if it could instanta- 
neously follow the modulation of the time-dependent 
parameters. Corrections due to retardation effects are 
encoded in p[ k> °\ and are governed by a competition 
between the time scales of the time-dependent driving 
and the response times contained in Wt- Since the ker- 
nel Wt depends parametrically on time (as emphasized 
by the subscript t), also the solutions of Eq.©, p[ k \ 
acquire a parametric dependence on t. 

We stress that the validity of Eq.© is restricted to 
the regime of weak coupling to the leads HT <IC fceT and 
slow driving. More precisely the latter condition requires 
that the frequency and the amplitude of the modulation 
satisfy the relation 5x& <C Tk^T, where Sx represents 
the amplitude of any of the modulation parameters. In 
this limit, the expansion in powers of for the occupa- 
tion probability can be truncated to k = 2, as it will be 
discussed throughout this paper. 



C. Charge and heat currents 

We are interested in the charge and heat currents flow- 
ing in the leads in the presence of an external time- 
dependent driving. For definiteness, we take the sign 
convention that in each lead the particle (heat) current 
is positive when flowing towards the double dot. With 
this convention, the charge current in lead a can be writ- 
ten as 



I a (t) = ejTr {N a p(t)} 



(6) 



where N a = J2ka c Lka Caka ' ^ s ^ ne occupation-number op- 
erator in lead a, and p(t) is the density matrix of the total 
system (double dot and leads) and e > is absolute value 
of the electron charge. Similarly the heat current in lead 
a is given by 



H a — p a N c 



Pit)} 



(7) 



and it represents the (negative) rate of change of the en- 
ergy in lead a measured with respect to the local chemical 
potential. 

Performing an expansion in the driving frequency 
similar to the one carried out for p(t), these currents can 



be expressed as a series of contributions of order tt k , 



k>0 



fc>0 



In the limit of weak coupling to the leads (fiT <C k-gT), 
to lowest order in T these are given b y 53 i 54 



J< fc )(i)=eJ>f\ 



(9a) 
(9b) 



where e = (1, 1, 1, 1, 1, 1, 1, 1, 1) and the kernels X" and 
J'l take respectively into account the charge and the 
heat that flow from lead a into the double dot 



-t kj 



-e(n, - rij) [W a ,iL 



uria - m - e s ) - ti^m - nj )} [w a A 



(9c) 
(9d) 



Here, and n-i are the energy and the number of elec- 
trons in the double dot in state \i), respectively. The 
matrix elements of the kernel, [W a ,t]ijj represent the 
probability per unit time that a tunneling event from/to 
lead a induces the transition \j) — > \i) (i ^ j), to lowest 
order in T. 

The zeroth-order contributions Ia^ (t) and Ja (t) rep- 
resent the steady-state charge and heat currents that 
would flow in the leads in a stationary situation with the 
time-dependent parameters "frozen" at time t. They are 
non zero only if the system is brought out of equilibrium 
by applying a bias voltage or a temperature gradient. In 
the following we will also refer to them as "instantaneous 
currents" , because they would be the only contributions 
to the currents if the system could follow instantaneously 
the time-dependent driving. 

Vice versa, terms with k > describe the additional 
contribution to the currents due to the delayed response 
of the system to the time-dependent modulation. From 
Eqs. © and the sum rule X)jfWt]ij = 0, it follows di- 
rectly that they satisfy the identities 



t\t)+l£\t) = -eUn ) r i \ 



dt 



(10) 



where (h)?> = E.n^J and /(*>(*) = |(jW(t) 
r(*0 



(fc)/ 



Jl (t)) are the corrections of order k to the occupation 
of the double-dot and to the current flowing through it. 

The first equation represents essentially the charge 
continuity equation, and it ensures the conservation of 
the charge at every order of the frequency expansion^ 
A direct consequence of Eq. (fTUf is that the time-averaged 
electric current is uniform = —1^ , where the bar 
indicates the average over one driving period 



2^ 



5 



The second equation expresses instead the first princi- 
ple of thermodynamics, which relates the total heat flow- 
ing into a system to the increase of internal energy of the 
system itself and the work done by the external sources. 
This becomes more evident by rewriting the first term 
on the right hand side of Ea. (fTTj) in terms of the inter- 



nal energy of the double dot {E)t = J2i EipfJ and the 
power delivered by the ac-sources applied to the gate 



Ik) 



.31 



-p(*0 _ 

/ ac — 

comes 



E a =L,R dt 



Ak) 



(fc-1) 



In this way Ea. (fTTj) bc- 



(fc-i) 



dt 



(ft) 

dc ' 



(12) 



where VfJ = VI^ is the power delivered by the external 

dc-source according to Joule's law. Note that V^. is 
positive if the current flows in the direction set by the 
bias, and negative otherwise (sec discussion in Scc. lIII CI) . 
Integrating such an expression over one cycle gives 



(ft) 



4 fc) = -(^ fe) 



(13) 



which allows to express the work per cycle done by the 
ac-sources in terms of the heat currents flowing in the 
leads. 



1.56 



III. CHARGE AND HEAT PUMPING 

The time-dependent modulation of the voltages ap- 
plied to the local gates of the quantum dot allows to 
pump charge and heat through the system. In the fol- 
lowing, we will discuss the transport features of such a 
double-dot pump and how it can be understood as a 
nanoscale engine. In all the calculations presented be- 
low, the pumping cycle is parametrized in terms of the 
mean energy E(t) = E + 5ECOs(ttt + (f>) and the detuning 
e{t) — e + 5 e cos(fif). 



A. Pure adiabatic pumping 

To understand the features in the heat and charge cur- 
rents due to time-dependent driving, we start by consid- 
ering the case of pure adiabatic pumping, meaning that 
the two leads are at equilibrium (Tj, = Tr = T and 
V = 0) and transport is only due to the slow modulation 
of the levels of the two dots (f2 — s- 0). In this case, the 
relevant contributions to the currents are the ones to first 
order in the driving frequency, Ia\t) and Ja\t). 

If the leads are in equilibrium, the charge current 
pumped through the double-dot takes the simple formal 



77— a, b 



with p^' 
where pf^ = e 



E<7=t4^° ' aild 



= T v {0) 
. It can be divided 
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FIG. 2. (Color online) Contributions to first order in Q to the 
average charge (a) and heat (b) currents, plotted as a function 
of the mean energy E. The pumping cycle is defined by: 
8 e = 2Se, £ — 0, (f> = 7r/2 and SI = T/200, and it corresponds 
to a circular orbit centered around zero detuning. In both 
panels: T L = T R = T, V = and U = 20k B T, t c = Wk B T, 

r L = r R = r/2, hr = fc B r/4. 



in two contributions /i 1 ' = Iyl 



each of which con- 
tains only transitions that involve tunneling in or out of 
the bonding or the antibonding state, respectively. These 
two terms always contribute to the current with oppo- 
site signal and, while the first is dominant around the 
triple point located at (e^e^) « (0,0), the second one 
is largest around (€l,€r) ~ (~U, —U). This leads to the 
sign change of the pumped charge as a function of the 
mean energy E shown in Fig.[2ri. Here, we plot the time- 
averaged charge current flowing into the right lead for 
various amplitudes of the ac-modulation applied to the 
gates that forces the state of the system to follow orbits 
in parameter space similar to those of Fig. [TJ;. When the 
amplitude is large enough so that the pumping cycle fully 
encircles a triple point, see e.g. the red orbit in Fig. [TJ;, 
electrons can be transferred one by one through the dou- 
ble dot, generating a quantized dc-current / = ±eQ/2ir. 

One of the key ingredients of such a quantized pump- 
ing regime is the alternate decoupling from the leads££ 
whenever one of the two dots comes in resonance with its 
neighboring lead, the other one is strongly off- resonant, 
so that particle exchange occurs only with one of the 
leads at a time. In Eq. (fH)) , the "alternate decoupling" is 
encoded in the time-dependent prefactors T a ^ v /T v and it 
can be achieved only if the detuning is much larger than 
the inter-dot coupling e ^S> t c at the time when the level of 
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one dot crosses the Fermi energy of the neighboring lead, 
see Eq.([3]). This in turn requires the modulation ampli- 
tude to be larger than t c . The second ingredient of quan- 
tized charge pumping is the strong inter-dot Coulomb 
interaction U, which separates the two triple points well 
apart and which forbids, for example, the double occupa- 
tion of the double dot along an orbit that encircles only 
the triple point around (e^,, cr) ss (0, 0). This permits to 
transfer charges sequentially from one lead to the other 
through the two dots according to the direction set by 
the pumping cycle. Vice versa, no directional transfer is 
possible if the orbit encircles both triple points, so that 
the maximal width of the plateaux in the pumped charge 
current is of the order U + A, further reduced by tem- 
perature smearing. 

The heat current in the leads has in general a rather 
complicated analytic expression, but in the limit in which 
the occupation of the anti-bonding state can be neglected 
(i.e. t c 3> k&T), it is well approximated by 

Ja £b r b dt Po + ^ + u > r a dt Pd ■ (15j 

This equation indicates that a change in the probability 
for the double dot to be empty or doubly occupied results 
in a heat current in the leads that is directly proportional 
to the energy involved in the transition responsible for 
the change. These two contributions arc weighted by dif- 
ferent time-dependent prefactors, r Qj b/rb and r Qja /r a , 
which results in the fact that, as for the charge current, 
changes in p^ and p^ contribute with different signs to 
the heat current. 

Interestingly, in the regime where the charge current 
is quantized, also the heat current shows well defined 
plateaux with height i^fesT In 2, see Fig. [5b. The emer- 
gence of these plateaux can be well understood. In fact, 
in the regime of quantized charge pumping, J„' can be 
directly related to the entropy difference AS'i" ' between 
two charge configurations that differ only by one electron 
in dot a: 

^ J« = k B T f^dt^lls^ = k B TAS^ . (16) 

" JO l rj at 

Here, = — Y^iPi m Pi is the Shannon entropy of 
the double-dot and 77 = b (a) for an orbit fully encir- 
cling the triple point around (£l,£r) ~ (0,0) (around 
(cl ) £r) ~ (— U, — U)). Along this orbit, the time- 
dependent prcfactor T a ^ ri /T rj equals one close to the res- 
onance with lead a and zero otherwise, so that the in- 
tegral results in the difference in entropy of the double 
dot before and after an electron has tunneled through 
the a barrier, i.e. ASi ' = ±ln2. The plus or minus 
sign corresponds to an electron tunneling in or out of the 
double-dot, according to the direction set by the pump- 
ing cycle. 

Note that fee Tin 2 is the minimal amount of energy 
required to erase one bit of information, according to 




-100 -50 50 



E [k B T] 

FIG. 3. Contributions to first order in fl to the average charge 
and heat currents for the case of a fully spin-polarized system. 
The pumping cycle is defined by: 8b = 40/cbT, 8 e = 2Se, 
s = 0, 4> = tt/2 and Q = T/200. Other parameters: T L = 
T R = T, V = and U = 20fc B T, t c = 10k B T, T L = r R = T/2, 
fir = fc B T/4. 

Landauer's principle^ - — In this case the bit of informa- 
tion is encoded in the spin of the electron in the bonding 
state, the erasing procedure corresponds to raising the 
level of one dot and allowing the electron to tunnel out. 
The energy /cBTln2 is provided by the external ac-ficlds, 
and results in heat that flows into the lead the electron 
has tunneled to. The situation is specular for the re- 
verse process, so that at the end of the cycle the energy 
/cbT In 2 has been transported from one lead to the other. 

While the appearance of plateaux in the charge current 
is directly related to the quantization of charge, plateaux 
in the heat current reflect specific degeneracies occurring 
in the system and are therefore tunable e.g. by an ex- 
ternal magnetic field. This is shown in Fig. [21 where 
the plateaux in the heat current are fully suppressed 
by a strong magnetic field that spin-polarize the system 
B ^> U-qT . Vice versa those in the charge current are 
essentially unaffected. In this case, charge transport is 
not accompanied by heat transfer from one lead to the 
other. From the information-theory point of view, one 
could say that the regular flow of spinlcss electrons pro- 
duced by the double-dot pump in a high magnetic field 
does not carry any useful information and, therefore, it 
is not accompanied by the heat transfer which otherwise 
would be required by Landauer's principle. 

Another feature that distinguishes heat from charge 
pumping is the appearance of peaked shaped features 
at the borders of the plateaux. These peaks approach 
the values ±^/cBTln3 and ±|^fcBTln| and are asso- 
ciated to orbits that touch one of the triple points, see 
e.g. the black orbit in Fig. [TJ;. In this case particles are 
also exchanged with the leads in situations in which two 
different charge states (e.g. singly and doubly occupied) 
are degenerate, and this results in differences of the en- 
tropy between initial and final state of a certain loading 
and unloading process that equal ± In 3 or ± In | . This 
change in entropy does not translate entirely in the heat 
exchanged with one lead because close to a triple point 
the energy dependence of the weight factor r a) ^/r^ in 
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Eg. (fTB|) . cannot be neglected or, equivalently, the "de- 
coupling approximation" breaks down and each dot ex- 
changes heat with both leads simultaneously^ 

Finally, we notice that if the leads are in equilibrium 
Eq.(JTTj) reduces to 

jl 1 \t) + 4\t) = k B T±sW. (17) 

This expression is equivalent to the Clausius equality of 
equilibrium thermodynamics and it indicates that to first 
order in fi the dynamics of the system is reversible if the 
leads are in equilibrium. A direct consequence of Ea. (fT7)) 
is that on average, during a cycle, heat is transported 
from one lead to the other — — • 

B. Pumping against a gradient 

So far we consider only the case of pure pumping, in 
which the leads are in equilibrium and transport is solely 
due to periodic modulation of dots' levels. However, an 
important feature of the double-dot pump is the possi- 
bility of achieving quantized charge pumping even in the 
presence of a finite bias voltage V. From the experimen- 
tal point of view, pumping against a bias represents a 
clear proof that the measured current is set by the cho- 
sen pumping cycled Moreover, it also allows thinking 
of applications of the double-dot pump e.g. as "battery 
charger" transferring electrons from a lower to a higher 
chemical potential. This will be discussed in the follow- 
ing Section UTTCl 

Pumping against a voltage bias requires an orbit in 
parameter space that minimizes the contributions of the 
instantaneous currents I a and J a , which - flowing in 
the direction set by the gradient irrespectively of the ori- 
entation of the pumping cycle - play the role of leakage 
currents. This can be achieved by choosing a pumping 
cycle that fully encircles a single triple point. However, 
it has to be taken into account that the triple points, i.e. 
the regions in parameter spaces where no charge configu- 
ration is stable, get broadened by the gradient itself and 
will eventually merge into a single one for very large V. 
This poses an upper limit to the maximal voltage against 
which it is possible to pump, see Sec. IIVBI Nevertheless, 
as long as the triple point regions are well separate from 
each other and can be encircled by a controlled modula- 
tion of the gates, it remains possible to pump one electron 
per cycle through the double dot. 

The situation is similar for what concerns heat pump- 
ing when the leads have different temperatures: as long as 
it is possible to choose a pumping cycle that fully encloses 
a triple point, the heat currents in each lead exhibits well 
defined plateaux of height ■^kBT a \n2, see Fig. |4l In 
other words, even if the system is globally brought out of 
equilibrium, the heat exchanged singularly with each lead 
obeys Clausius relation J a = /sbIqAs! ^. This is 
again a consequence of the "alternate decoupling" occur- 
ring along an orbit fully encircling a triple point, which 
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FIG. 4. (Color online) Contributions to first order in Q to the 
average heat current in the left (dashed lines) and the right 
lead (full lines), with and without a temperature gradient. 
Here we assumed T L = T+ AT and T R = T, with AT = (red 
lines) or AT = 0.4T (green lines). The thin dot-dashed lines 
correspond to <;b(T+AT) In 2Q./2n. The heat current in each 
lead shows plateaux proportional to the local temperature 
T a . When j£ < and J^ 1 ' > 0, heat is pumped form the 
cold to the hot reservoir . The pumping cycle is defined by: 
5 E = 40fc B T, 5 e = 25 E , e = 0, <j> = vr/2 and 9, = T/200. 
Other parameters: U = 20k B T, t c = 10fc B T, T L = Ir = T/2, 

hT = fc B r/4. 



ensures that the double dot is temporarily coupled only 
to one lead at a time. In the considered limit of slow 
driving O — t 0, the double-dot has time to equilibrate 
with the lead it is coupled to. As a consequence, even 
in the presence of a temperature gradient, processes that 
change the total occupation of the double dot represent 
isothermal transitions between equilibrium states. 

In the limit of slow driving, a pumping cycle that 
fully encloses a triple point can then be regarded as 
a nanoscale analog of the Carnot cycle in which two 
isothermal transitions are connected by two adiabatic 
ones (where now adiabatic means that the entropy of the 
system remains unchanged). The isothermal processes 
are those that change the total occupation of the double 
dot, and the adiabatic ones those that accompany the 
crossing of the two dots' levels, which, occurring when 
both dots are far away from the Fermi energy of the leads, 
represent reversible processes. 

Finally we notice that the total average heat-current 
flowing into the system j£ + Jr is directly related to 
the excess entropy production due to the time dependent 
driving^ As emphasized in Refl62l. this is a geometric 
quantity that depends only on the trajectory in param- 
eter space. The geometric nature of the excess entropy 
production - and more in general of the pumped charge 
and heat to first order in frequency - is a direct con- 
sequence of the relation between and Ja and the 
time derivative of the instantaneous occupation proba- 
bility , see Appendix. IAW 
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C. Operation of a double-dot device 



2. Heat pump 



The possibility of pumping charge and/or heat against 
a gradient shows that a driven double-dot device can 
have different applications as a nanoscale engine. In the 
following we discuss in detail three different operating 
modes: the charge pump moving electrons against a bias 
voltage, the heat pump transferring heat from a cold to 
a hot reservoir, and the thermal engine producing work 
extracting energy from a hot reservoir and dumping part 
of it in a cold one. To quantify their performance, we 
introduce efficiencies in analogy to the classical counter- 
parts. 

In this section we concentrate on the ideal working 
regime where the main contributions to the pumped 
currents are of first order in the driving frequency, i.e. 
I a fa Ia^ and J a « J„ . Corrections due to finite 
operation-time and leakage currents will be analyzed in 
Section HVl 



1. Charge pump 

When the double-dot is operated in a way as to transfer 
electrons from a lower to a higher chemical potential, it 
can be regarded as a "battery charger" doing work on 
the dc-source. We characterize the performance of such 
an engine by its energy conversion efficiency 



As discussed in Sec. IIIIB1 a double dot operating be- 
tween two reservoirs with different temperatures along 
an orbit that fully encloses a triple point, represents an 
analog of the Carnot engine. Depending on the direction 
of the cycle, the device acts either as a refrigerator or as 
a heat engine. We are interested here in the situation 
where no voltage bias is present. 

The double dot is operated as a refrigerator when heat 
is transferred from a cold to a hot reservoir during the 
pumping cycle. The efficiency of a heat pump is in gen- 
eral characterized by the coefficient of performance 



COP 



J cold 



(20) 



where J co id is the average heat flow out of the cold reser- 
voir and P ac the power delivered by the ac-field during 
one cycle. If no bias is applied, the latter is directly re- 
lated to the total heat flowing out of the double-dot P ac = 
— J2 a 4i see Eq. ([T5)l . Along an orbit that fully encircles 
a triple point, the heat current in each lead is directly 
proportional to the temperature of that lead. In this case 
Xoid = £fc B T co i d ln2 and P ac = £k B {T hot - T cold )ln2 
and the cooling efficiency of the double-dot pump reaches 
the Carnot limit COP Ca rnot = T cool /(T hot - Tcooi)) con- 
firming the analogy between the "ideal" pumping cycle 
and the Carnot cycle. 



Vch = 



-IV 
V : 

' ac 



(18) 



3. Heat engine 



where —IV is the useful work per unit time done by the 
pump on the dc-voltage source, and P ac is the average 
power delivered by the ac-fields applied to the gate elec- 
trodes. The sign convention used here for the current is 
that / is positive when it flows in the direction set by the 
dc source, so that —TV > when the driven quantum 
dot actually works as a charge pump moving electrons 
from a lower to a higher chemical potential. Only in this 
case it make sense to speak about useful work done by 
the pump and to characterize its efficiency in terms of 
Eq.CEH). 

The power delivered by the ac-fields V ac can be quan- 
tified in terms of the heat currents flowing into the leads 
according to Eq. ([lU)) . i.e. P ac = — P dc — J2 a Ja, where 
Vdc = IV is the work per unit time done by the dc- 
source. The efficiency of the double-dot pump as a bat- 
tery charger is then given by 



IV 



Jl + Jr + IV 



(19) 



If the leads have the same temperature, the total heat 
current is zero along an orbit that fully enclose a triple 
point and the efficiency takes its maximal value ?7 c h = 1. 
In this case, the work done by the source is fully used to 
transfer charge against a potential. 



The double-dot pump acts as a heat engine producing 
work on the external fields, when the pumping cycle is 
such that heat is extracted from the hot reservoir and 
released into the cold one. For the case of Fig. |4l this 



corresponds to the regimes where > and J^> < 0. 
The performance of a heat engine is characterized by the 
efficiency coefficient 



V 



-P a 



(21) 



where — V ac is now the work per unit time done by the 
system on the ac fields and Jhot is the heat current ab- 
sorbed from the hot reservoir. 

Similarly to the case of the refrigerator, along an orbit 
that fully encircles a triple point it is —P ac — ^-feBC^kot - 
7coid) m 2 an d Jhot = ^-fcB?hot hi2, and the efficiency 
of the double dot pump reaches again the Carnot limit 

^Carnot = 1 _ 7c ld/?hot- 



IV. LIMITATIONS OF A REALISTIC PUMP 

In the previous section we have assumed that the 
charge and heat currents flowing in the leads can be ap- 
proximated solely by the contributions to first order in 
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FIG. 5. (Color online) Contributions to second order in SI 
to the average charge (a) and heat (b) currents, plotted as a 
function of the mean energy E. The pumping cycle is defined 
by: 5 £ = 26 E , £ = 0, cj> = it/2 and SI = T/200. The thin 
lines in panel (b) represent the approximate expression Ja fH 
In both panels: T L = T R = T, V = and U = 

20k B T,\ = io/cbT, r L = r R = 172, nr = fc B r/4. 



frequency, I a w 1^ and J a 



J (1) 



This correspond to 



considering the case of extremely slow driving fl — > and 
to implicitly assuming that the instantaneous currents 
Ia^ and Ja^ can be neglected along an orbit that fully 
encircles a triple point. In this ideal case, the efficiencies 
of the double-dot based engines discussed in Sec. MI CI 
reach their maximum theoretical values. In this section 
we investigate instead the limitations to the performance 
of a realistic double-dot pump that are due to a small 
but finite driving frequency and to leakage currents. 



A. Corrections due to finite operation-time 

To take into account the effects of a finite driving fre- 
quency, we evaluate the contribution to second order in 

(2) (2) 

O to the charge and heat currents I a , Ja ■ For simplic- 
ity, we start again by first discussing the case where the 
leads are at equilibrium and postpone the case of pump- 
ing against a gradient to Sec. IIVBI Moreover, we focus 
on the regime of large modulation amplitudes, which is 
the most interesting for applications. The opposite case 
of small amplitudes pumping is described in detail in the 
Appendix. [B] 

In Fig. [5] we plot the contributions to second order in 
SI to the average charge and heat currents I R , J R , for 
the same choice of parameters considered in Fig. [21 Com- 
paring these figures it is apparent that, even in the case 
of relatively slow driving SI = r/200 considered here, the 
contribution to second order in SI to the heat current rep- 
resents a significant amount of the total heat exchanged 
with a lead per cycle. Vice versa, the contributions to sec- 
ond order in to the charge current can be completely 



neglected. 

This behavior is easily understood considering again an 
orbit in parameter space that fully encloses a triple point. 
As discussed before, along such an orbit the exchange of 
particles and heat with the leads involves only one dot 
and its neighboring lead at a time, while the other lead 
is temporarily decoupled due to the large inter-dot de- 
tuning. Particle exchange occurs when the level of a dot 
crosses the Fermi energy of the neighboring lead. How- 
ever, if the level moves with a finite velocity, electrons 
can tunnel out (in) of the dot not only in resonance, 
but also at a somewhat higher (lower) energy, resulting 
the emission of in a hot electron (hole) into the leads. 
Both processes contribute to increasing the energy of the 
lead and give rise to a negative heat current (meaning 
that it is flowing from the dot into the leads). In other 

—(2) 

words, J a represents the irreversible heat production 
that accompanies the action of the ac-fields applied to 
the double dot. The maximal excess energy that can be 
deposited into the leads is determined by the amplitude 
of the modulation 8e, while the probability that strongly 
off-resonant transition occurs is given by the ratio S7/r, 
so that the heat current to second order in SI scales as 
\Ja 2 \ ~ fe^ 2 /r. A more accurate estimate leads to 



Ji 2 



3 v a SI 



(22) 



where v a is the speed with which the level of dot a crosses 
the Fermi energy of lead a. For a circular orbit fully 
encircling the triple point around (cl,£_r) = (0,0) it is 
v a = Q\/25e\/1 — e^/2S E ^ For an orbit around the 
other triple point, e a should be replaced by (e a + U) 
in the expression for v a . Deviations from the behavior 
predicted by Eq. (|2"2"j) occur for orbits that cross a triple 
point and are associated to an increased heat production 
close to the turning points of the ac-modulation. The 
large peak around E = —U/2 for large amplitudes corre- 
sponds to orbits that enclose both triple points. In this 
case a dot level crosses the Fermi energy of the neighbor- 
ing leads twice during each cycle, leading to (almost) a 

(2) 

doubling of the heat current J a . 

It is important to stress that the contributions to 
second order in S7 to the heat current are negative in 
both leads (with = in the case of zero aver- 
age detuning e = 0)^ meaning that heat is deposited 
in the leads in every pumping cycle. This pose some 
serious constrains to the possibility of using a double- 
dot as a heat pump to extract heat from a reservoir, as 
it requires to reduce significantly the driving frequency 
S7 <C hi2TkBT/SE in order to minimize heating effects. 

While the emission of a hot hole or a hot electron con- 
tribute with the same sign to the heat current into the 
leads, the two processes compensate each other in the 
average charge current, so that I a remains in general 
pretty small as long as SI <C T. In particular, it can 

be shown that I a = for orbits with average zero de- 
tuning, thanks to the symmetry of the pumping cycle, 
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FIG. 6. (Color online) Average charge current as a function of 
the mean dot-level energy E for different values of the average 
detuning e. The thick lines represent the total charge current 
In = 1^ + while the thin-dashed ones correspond to 
the contribution to first order in frequency 1^ alone. The 
pumping cycle is defined by: Se = 40/cbT, S e — 2Se, s — 0, 
4> = tt/2 and Q = 17200. Other parameters: T L = T R = T, 
V = and U = 20fc B T, t c = Wk B T, I\ = Ir = 172, KT = 
ksT/4. Curves are offset by |ef2/27r for clarity. 

sec Fig. Even for orbits with e ^ 0, the corrections 

-(2) 

due to la to the charge current are negligible as long 
as the pumping cycle fully encircles a triple point, and 
become sizable only at the edges of the current plateaux, 
see Fig. [5] 

Finally, we would like to comment briefly on the consis- 
tency of our expansion in SI. As explained in Scc lIIBl our 
calculations are based on the expansion p(t) = J2k Pt > 
where the contributions pj; scales as (d,/T) k giving a 
convergent series in the limit ~ T. The scaling of the 
various contributions to the average charge and heat cur- 
rents is less obvious, as shown by J a 2 > J a in certain 
frequency regimes. This is because the physical processes 
responsible for the heating induced by the ac-driving 
scale as Jl 2 and therefore start contributing to the heat 
current only from second order in frequency, in terms of 
the expansion we use. However, it can be shown rigor- 
ously that all the contributions J a k ~ 2 ^ are suppressed at 
least by a factor f2/r with respect to J a , so that it is 
meaningful to truncate the expansion in Ea. (|9b[) to sec- 
ond order only. 



B. Leakage currents & efficiency of a realistic pump 

The leakage currents I a °^ and J a ^ pose another im- 
portant limitation to the operation of the double-dot as 
a heat or charge pump. In fact, it is possible to trans- 
fer electrons against a bias only as long as the pumped 
current I^p + I a 2 \ i.e. the current resulting from the 
time-dependent modulation of the dots' levels, prevails 

over the instantaneous one Ia ^ ■ Similarly, to pump heat 
from a cold to a hot reservoir it is necessary that the 



magnitude of J a ^ + Ja 2 " 1 is larger than the one of J a 0) . 
The contributions from I a ^ and J a °^ can be strongly sup- 
pressed by choosing an orbit along which the levels of the 
two dots cross each other only well outside of the bias or 
temperature window. However, for a given orbit, the 
maximal amount of charge or heat that can be pumped 
in or out of a reservoir in one cycle is fixed, while the 
leakage currents increase with V and AT. This poses 
an upper limit to the maximum voltage or temperature 
gradient against which it is possible to pump. Moreover, 
the contribution of I a ^ and J a ^ to the total charge and 
heat extracted from a lead during a pumping cycle grows 
as Q^ 1 when the driving frequency is reduced. As a con- 
sequence, the maximum voltage or temperature gradient 
that the double-dot pump can overcome depends both 
on the chosen pumping orbit and on the speed at which 
it is traversed. 

The competition between pumped and leakage currents 
can be seen e.g. in Fig. [7k,, for the case of a chargc-pump 
working against a bias voltage. Here we show the to- 
tal charge current I = 1^ + I^ 1 ' + 1^ flowing through 
the dot as a function of the applied bias for different 
amplitudes of the driving cycle and different driving fre- 
quencies. In each case it is 1^ « 0. According to our 
convention on the sign of J, the pump moves electrons 
against the direction set by the bias as long as 7 < 0. 
The maximal voltage V max that the pump can sustain is 
then defined by the condition 7 = 0. For larger voltages, 
the pumped current is overrun by the one due to the bias. 
As expected, Vm ax is sensibly reduced if the radius of the 
orbit is shrunk, or if the pumping frequency is lowered. 

As long as / < 0, the pump performs a positive work 
per pumping-pcriod, by transferring electrons from a 
lower to a higher chemical potential, sec Fig. [7)d. The 
efficiency of such a pump can be quantified by the energy 
conversion coefficient r] c h, Ea. (|19p . where now all contri- 
butions up to second order in to the charge and heat 
currents have to be taken into account, see Fig. [7fc. At 
large voltages, the major limitations to the performance 
of the pump come from the leakage current 1^ , which 
poses an upper limit to the maximal power that can be 
delivered by the pump on the dc-source. Vice versa, at 
low bias the efficiency r) c h is mostly affected by the heat 

-(2) ~(2) 

production term + , which represents the mini- 
mal power that the ac-sources have to provide in order to 
drive the system along the considered orbit. Since this is 
finite even in the limit of no applied bias, the efficiency 
of the double-dot charge pump vanishes for V —> 0. In- 
terestingly, despite the detrimental effects due to heat 
production and leakage currents, the double-dot charge 
pump can reach efficiencies up to 70% of the ideal value, 
when operated against a finite voltage along a sufficiently 
large pumping cycle, see Fig.. [7jp. 

The limitations due to leakage currents and finite 
operation-time are more severe when the double-dot is 
employed as a refrigerator to pump heat from a cold to a 
hot reservoir. In fact, while the contributions to second 
order in frequency are mostly negligible for what concerns 
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FIG. 7. (Color online) (a) Average charge current flowing 
through the double-dot in the presence of a bias voltage and 
of time-dependent driving. The current is plotted as a func- 
tion of the applied bias V for three different sets of pump- 
ing parameters: Curves (1) and (2) correspond to the orbit 
described by 5 E = 40fc B T, E = 30fc B T, e = 0, (f> = tt/2, 
traversed with frequency fl — T/200 and = T/400, re- 
spectively. Curve (3) correspond to a smaller orbit with 
S E = 20fc B T, E = 15fc B T, e = 0, <f> = tv/2 and Q = T/200. 
Charge is pumped against the bias as long as I < 0. (b) 
Work per unit time done on the dc-source and (c) efficiency of 
the double-dot pump for the same pumping cycles considered 
in panel (a). The efficiency is defined according to Eq. (|19|) . 
In all panels: U = 20fc B T, t c = Wk B T, F L = T R = T/2, 
hV = fc B T/4 and T L = T R = T. 



the charge current, they represent a significant amount 
of the heat flowing into the leads. Moreover, being al- 
ways injected from the double-dot into the leads, they 
reduce the possibility of extracting heat from a reservoir. 

~(2) 

Suppressing Jq requires to lower considerably the driv- 
ing frequency, but this in turn makes the system more 
sensitive to leakage currents. Both effects reduce signif- 
icantly the efficiency of the double-dot pump with re- 
spect to the Carnot limit. This can be seen in Fig. [5^,, 
where we plot the coefficient of performance of a double- 
dot heat pump, plotted as a function of the temperature 
gradient between the leads for the same pumping cycles 
considered in Fig. [7J The Carnot efficiency is also plotted 
for comparison. The effects due to finite-time operation 
are particularly evident in the limit AT — > 0, where an 
ideal pump would be "infinitely efficient" , since it re- 
quire no power to transfer adiabatically heat from one 
lead to the other when these have the same temperature. 
Vice versa, due to the heating of the leads which accom- 
panies the modulation of the levels of the double-dot, 
the power needed to operate a pump at finite speed re- 
mains finite, as well as its efficiency. Moreover, an ideal 
heat pump would be able to work against an arbitrary 
temperature gradient, but this is not the case of a real 
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FIG. 8. (Color online) (a) Coefficient of performance of the 
double-dot pump as a heat pump, as defined in Eq. (|20[l , In 
this case we choose TL = Tr + AT and a pumping cycle such 
that heat is transfered from right to left, i.e. against the tem- 
perature gradient. The curves labeled as (l)-(3) correspond 
to the pumping cycles discussed in Fig. [7] COP of an ideal 
Carnot refrigerator is also plotted for comparison, (b) Effi- 
ciency of the double-dot pump as a heat engine, as defined in 
Ea.lpl)). Also in this case TL = Tr + AT, but the pumping 
cycle is now such that heat is transfered from left to right, i.e. 
according to the temperature gradient. The curves labeled as 
(l)-(3) corresponds again to the orbits discussed in Fig. [7J but 
now traversed in the opposite direction, i.e. <j> — —n/2. The 
efficiency of an ideal Carnot engine is plotted for comparison. 
In all panels: U = 20fc B T, t c = 10fc B T, T L = Tr = T/2, 
fir = fc B T*/4 and /il = /iR. 



double-dot pump, where the pumped current is overrun 
by leakage currents as soon as the temperature gradient 
AT becomes larger than a critical value. 

Heat production limits also the operation of a double- 
dot pump as a heat engine, see Fig.[5jD. In fact, while the 
maximal amount of work that the pump can perform ex- 
tracting energy from a hot reservoir and dumping it into 
a cold one is proportional to the temperature difference, 
the energy cost of moving along an orbit in parameter 
space - which is represented by the heat dissipated into 
the leads - is roughly independent of AT. This means 
that if the leads have similar temperatures, the work put 
into the system by the external ac-fields exceeds the one 
that can be extracted form the heat reservoirs, leading 
to a negative energy balance. In this case, the pump pro- 
duces no work and it make no sense to speak about its 
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efficiency in terms of Eq. ([TT)|) . Dissipation due to heat- 
ing can be reduced by driving the system at a lower fre- 
quency, but at the cost of increasing the effects of the 
leakage currents. These drag a substantial amount of 
heat out of the hot reservoir without contributing to the 
work performed by the system (since + J^ ' = 0), 
suppressing the efficiency of the pump for large AT. 

Finally we would like to mention that, apart for leak- 
age currents and heat production, other effects, such as 
coupling to phonons, cotunneling and missed inter-dots 
transitions, may affect the performance of the double-dot 
pump. However, the last two mechanisms are suppressed 
in the limit of weak coupling to the leads, strong inter-dot 
hybridization and slow driving frequencies considered in 
this paper, and coupling to phonons is also expected to 
be significantly reduced at the cryogenic temperatures at 
which pumping experiments are performed. 

V. CONCLUSIONS 

We investigated charge and heat transport in a driven 
double-dot device. We showed that in the regime of 
quantized charge pumping, i.e. when electrons are trans- 
ferred one by one through the system thanks to a con- 
trolled modulation of the energy levels of the two dots, 
the heat current exhibits well defined plateaux if the driv- 
ing frequency f2 is sufficiently small. The height of these 
plateaux is directly proportional to the temperature of 
the leads and reflects specific degeneracies of the double- 
dot states involved into transport, namely the spin de- 
generacy of the bonding state. In the limit of slow driving 
SI — >■ 0, the heat current through the double dot can then 
be controlled by an external magnetic field and be com- 
pletely suppressed by fully spin polarizing the system. 

While the quantization of the transported charge is 
rather robust with respect to an increase of the pumping 
frequency (as long as the latter remains much slower than 
the relaxation time of the system Q <C T), the plateaux in 
the heat current are strongly affected by heat production 
in the leads. This is described by the contributions to 
second order in the pumping frequency to the heat cur- 

(2) 

rents Ja , which account for the excess energy deposited 
into the leads as a consequence of the finite-time driving 
of the dots' levels. 

Heating effects are detrimental for the performance of 
double-dot pump as a nanoscale engines. In fact, while 
the driven double-dot can in principle be operated as 
a nanoscale "battery charger" moving electrons from a 
lower to a higher chemical potential, or as a heat pump 
exchanging energy with two reservoirs with different tem- 
peratures, its efficiency is limited by dissipative effects 
due to leakage currents and finite-frequency operation. 
We show that despite these effects, the efficiency of a 
double-dot charge pump performing work against a dc- 
source can reach of up to 70% of the ideal value. 

Our results are based on a real-time diagrammatic ap- 
proach in the regime of weak coupling to the leads, which 



allows to take systematically into account the effects of 
a small but finite driving frequency. It is interesting to 
observe how fundamental thermodynamic relations, see 
e.g. Eq. (fT2"j) . naturally emerge from transport calcula- 
tions and how the expansion in frequency allows to iden- 
tify reversible and irreversible transport processes. 

In this work we focused in particular on the semi-classic 
regime of strong inter-dot coupling and large modulation 
amplitudes, which is the most interesting for applications 
of the double dot as a quantized charge source or, as 
we discussed above, as a nanoscale engine. As future 
perspective, it will be interesting to extend our formalism 
to take into account coherences^ between the double-dot 
states and cotunneling effects^ to address regimes where 
quantum effects play a major role for the dynamics of the 
system. 
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Appendix A: Geometric properties of 1^ and J^ 1 ' 

It is well known that adiabatic pumping is of geomet- 
ric nature, meaning that, in the limit of slow modula- 
tions f2 — > 0, the charge pumped through a system de- 
pends only on the specific shape of the path sustained 
by the system's parameters but not on its detailed time 
evolution^S The same holds true for the heat ex- 
changed with the leads during one pumping cycle, as we 
discuss in this appendix. 

In terms of the expansion in frequency Eq. ([8]). the adi- 
abatic heat and charge currents correspond to the contri- 
butions to first order in the modulation frequency I a 
and J a . To emphasize the geometric nature of the 
charge and the heat pumped from a lead during one cycle, 
it is useful to express them in terms of auxiliary vector 
fields in the space of the parameters^ The key observa- 
tion is that la and Ja* 1 are directly related to the time 
derivative of the instantaneous occupation probabilities, 
i.e. 

^(t)=eX?W7 1 |pf ^fl a j/t\ (Al) 

and similarly for the heat current Ja\ Here, VV t is a 
pseudo-inverse of the evolution kernel Wt The second 
identity in Eq. (|Al[) corresponds to the definition of the 
response coefficients tp\, where £ is either I a or J a for 
the charge and heat current, respectively. 

The charge and the heat exchanged with one lead dur- 
ing a pumping cycle can then be written as a line integral 
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in parameter space 

r 2ir/n 



a) 



(A2) 



The vector e = J2 a e a^ a indicates the state of the sys- 
tem along the closed trajectory C in the parameter space 
spanned by eL = (1, 0) and = (0, 1). The vector field 



(A3) 



with V = ^2 a d £a e a , can be interpreted as a pscudo 
vector potential defined in the space of the driving pa- 
rameters. Its components describe the amount of charge 
(£ = la) or beat (£ = J a ) entering the leads due to 
a change in the driving parameters and are related to 
the emissivitiesi 68 i 69 The line integral on the right hand 
side of Ea. (|A2|) is independent of how fast the orbit C 
is traversed (provided that the adiabatic approximation 
holds), and therefore the pumped charge (or heat) does 
not depend on the driving frequency. 

Extending the two-dimensional parameter space to a 
third dimension and using Stoke's theorem, Eq. (|A2p can 
be written as the surface integral 



dS-B^e), 



(A4) 



where B^(e) = V x ^(e), S is the area encircled by C 
and dS the directed surface element in parameter space. 
The pumped "charge" can then be seen as the flux 
generated by the pseudo magnetic field B^. The advan- 
tage of this representation is that B^ anticipates the con- 
ditions for finite pumping without referring to the specific 
details of the modulation^ 

We stress however that the geometric interpretation 
of the pumped charge (heat) holds only to first order 
in the driving frequency and cannot be generalized to 
higher-order contributions, which depend sensitively on 
the details of how the pumping orbit is traversed. 



Appendix B: Small-amplitude pumping in first and 
second order in fl 



In order to get a better understanding of the charac- 
teristics of the pumped charge and heat, we present in 
this appendix some results in the regime of weak pump- 
ing. This means that the amplitudes of the parameter 
modulation are chosen to be small with respect to the ex- 
tension of the triple points. We furthermore take the am- 
plitudes of the time-dependent part of both energy levels 
= £l + <5el sin(Ot) and e R = l R + Se R sm(Qt + 4>) 
to be equal, Se^ = Se R = Se^- Finally we are inter- 
ested in the parameter cycles for which adiabatic charge 
and energy pumping is optimal and we therefore choose 
the phase difference between the parameters to be ir/2. 
These conditions result in modulation cycles of circular 
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FIG. 9. (Color online) Density plots of pumped charge and 
heat currents to first and second order in the pumping fre- 
quency, as a function of the energy level of the left and right 

dot. (a) Charge current to first order in Q,, 1^ [e£l ^ S ^yi ]i (b) 
Heat current to first order in SI, [Q %~pr ]> ( c ) Charge cur- 
rent to second order in f2, J^ 2 ' [eQ ^ S ^yi ], and (d) heat current 

= T, 

r/2, 



to second order in f2, J^ 2 ' [f2 ^ c j ]■ In all panels: TL = Tr = T 



V = and U 
HT = k B T/4. 



20fc B T, t c 



Wk B T, r L 



shape in the stability diagram. The contributions to the 
pumped charge in the small-amplitude regime in first or- 
der in the frequency are proportional to the area of the 
enclosed cycle 



2ir/f2 



e L (t)e R (t)dt = - 



2-n/n 



e L {t)e R (t)dt = irSe 



d- 



In order to further simplify our considerations we assume 
symmetric coupling, Tl = T R between the dots and their 
neighboring leads. Asymmetric couplings lead to quali- 
tatively similar results as the ones presented here, except 
for slight changes in the symmetry of the pumped charge 
and heat as a function of the working point (cl, £r). 

The results for the pumped charge and energy current 
in first order in the frequency are shown in Fig. [9£i,b. 
Both of these currents are pure transport properties, in 
the sense that we have = — T R as well as = 

—J R - Both charge and heat currents have contributions 
in the vicinity of the triple points in the stability diagram 
of the double-dot. The pumped charge changes sign be- 
tween the two triple points, which is due to the sign dif- 
ference in the dependence on the detuning of the effective 
coupling of the two hybrid states, T b / a , see Eq. 

The heat current has a different behavior: the contri- 
butions to the heat current have a node in the central 
regions of the triple points. This is due to the fact that 
the transported energy becomes zero when the quantum 
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dot level through which transport takes place is exactly 
at resonance. Furthermore we note that there is an asym- 
metry in the magnitude of the positive and the negative 
contribution. This is due to the fact that the maximum 
contribution to transport is shifted by a temperature- 
dependent factor, fcBTln2, with respect to the zero of 
the transported energy. 

Finally we note that in regime of small pumping- 
amplitudes, the average charge and heat currents to first 
order in f2 are directly proportional to the intensity of 
the pseudo-magnetic fields introduced in AppcndixlAl i.e. 
Ia oc Bi a and Ja^ oc Bj a , so that the plots of Fig. |9^,,b 
represents a color-scale map of Bi R and Bj R for the case 
V = and T L = T R . 

We now consider the second order in fi contribution to 
charge and heat transport. The pumped charge current 
1^ is - due to charge conservation - again a purely trans- 

— (2) — (2) 

ported quantity with 7^ = — 1^ . In contrast to the 
current in first order in 57, this contribution is symmetric 
as a function of the mean energy E around E = — y , and 
it is antisymmetric as a function of the detuning e around 
zero detuning. Another important difference is that it is 
due to single-parameter pumping, indeed we find 




Charge pumping due to a single time-dependent param- 
eter is not possible in the adiabatic regime. The different 
behavior of the second-order in frequency contribution 
can be understood by considering in detail a charge trans- 
fer process. Let us concentrate on the situation where the 
bonding level is close to resonance. The energy of the 
bonding level is decreased whenever the absolute value of 
the detuning is increased (the opposite effect is true for 
the antibonding state). A thereby induced loading of the 
bonding state takes place either mostly through the left 
or mostly through the right lead depending on the effec- 



tive coupling r L / R b . The fast change of the detuning 
has the consequence that the effective coupling "seen" 
by the tunneling particles is slightly delayed, leading to 
a slightly increased (decreased) coupling to the left lead 
with respect to the equilibrium situation when the de- 
tuning is increased (decreased), and vice versa for the 
coupling to the right lead (the opposite effect is true for 
the antibonding state). This allows for the appearance 
of single-parameter pumping, as reflected in the current 
formula, Eq. (|B1|) . However, since this delay in the effec- 
tive couplings is equal for a detuning of equal magnitude 
and opposite sign, charge pumping takes place only if the 
average detuning e is different from zero, see Fig. [§b. 

Finally, we want to address the heat current in sec- 
ond order in the frequency, which is due to the finite- 
frequency operation. This heat current is generated by 
the time-dependent fields and therefore has the same neg- 
ative sign in both leads. In other words, due to the time- 
dependent modulation, heat is flowing from the double- 
dot system into the leads. In Fig. [5U, the heat current 
into the right lead is shown. It shows large contributions 
whenever the level through which transport takes place 
(namely the bonding level at the transition H 1 and 
the antibonding level at the transition 1 <R- 2) is in res- 
onance with the right lead. Since the heat current in 
second order in the frequency is not a transported quan- 
tity in this regime, contributions also far away from the 
triple points are visible, as long as one of the levels is 
close to resonance. The analogous situation, with an in- 
verted behavior with respect to the detuning, is observed 
for the heat current into the left lead. 

This obviously different behavior of the heat current in 
second order in the frequency shows that here the second- 
order contribution is often larger than the first-order in 
frequency contribution. It is therefore of high relevance 
to consider contributions in second order in the frequency 
when studying heating effects. We however want to em- 
phasize that the rigorous expansion in the driving fre- 
quency remains well justified, since all higher contribu- 
tions in frequencies can be shown to be suppressed in 
powers of Q/T. 
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